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In the asymptotic limit Q 2 3> m 2 , the non-power corrections to the 
heavy flavour Wilson coefficients in deep-inelastic scattering are given in 
terms of massless Wilson coeffcients and massive operator matrix elements. 
We start extending the existing NLO calculation for these operator matrix 
elements by calculating the 0(e) terms of the two-loop expressions and 
having first investigations into the three-loop diagrams needed to 0(a;?). 

PACS numbers: PACS numbers come here 



1. Introduction 

In deep-inelastic scattering, the differential cross-section with respect to 
the Bj or ken- variable x and the virtuality of the photon Q 2 , can be expressed 
in terms of the unpolarized structure functions F%(x, Q 2 ) and Fl(x, Q 2 ), and 
the polarized structure functions gi(x,Q 2 ), g2(x,Q 2 ). For small values of 
x, the contributions of heavy charm product to F2{x,Q 2 ), F^x, (J 2 ), are 
of the order of 20 - 40 %, and therefore deserve and need a more detailed 
investigation. So far there exist NLO - 2-loop - heavy flavour corrections to 
F2 ,d (x, Q 2 ) in the whole kinematic range, calculated in a semi-analytic way 
in x-space [1]. A fast implementation for complex Mellin A-space was given 
in [2]. One observes that F^{x, Q 2 ) is very well described by an asymptotic 
result for F| c (x, (5 2 )|Q2 >m 2 for Q 2 >Wm 2 . For these higher values of Q 2 , 
one can calculate the heavy flavour Wilson coefficients, the perturbative 
part of the structure functions -^(x, Q 2 ) and Fl(x,Q 2 ), analytically, which 
has been done for F2(x,Q 2 ) to 2-loop order in [3,4] and for Fl(x,Q 2 ) to 
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3-loop order in [5] . First steps towards an asymptotic 3-loop calculation for 
F^^x, Q 2 ) have been done by the authors by calculating the first 0(e) terms 
of the 2-loop diagrams [6], contributing to 3-loop heavy flavour Wilson 
coefficients via renormalization. We report here on further steps towards a 
full 3-loop calculation. 



2. Heavy flavour Wilson coefficients in the limit Q 2 ^> m? 

On the twist— 2 level, the structure functions can be expressed as a convo- 
lution of perturbatively calculable Wilson coefficients and non-perturbative 
parton densities. We consider here the heavy flavour contributions to these 
Wilson coefficients, the heavy flavour Wilson coefficients. In the region 
Q 2 ^> m 2 , one can use the massive renormalization group equation to ob- 
tain all non-power corrections to these heavy flavour Wilson coefficients as 
convolutions of massless - light - Wilson coefficients Ck(Q 2 / ' n 2 ) and mas- 
sive operator matrix elements (OMEs) Aij{^ 2 /m 2 ) [1]. The light Wilson 
coefficients are known by now up to three loops [7] and carry all the pro- 
cess dependence. The operator matrix elements, on the other hand, are 
universal, process-independent objects, which are calculated as flavour de- 
composed operators in light-cone expansion between partonic states. Both 
objects have an expansion in a s , where one should note that the external 
legs of the diagrams contributiong to the OMEs are on-shell. 



3. Massive operator matrix elements 



In order to perform the 3-loop calculation of the OMEs, one has to 
first calculate the bare quantities and then to renormalize them, where they 
need to be mass- and charge-renormalized and contain ultraviolet (UV) 
and collinear divergences. The mass renormalization is done in the on- 
shell scheme [8,9], whereas the charge renormalization is done using the 
MS scheme. After mass- and charge-renormalization, the remaining UV- 
divergences are accounted for by operator renormalization via Z-factors, 
Zij, and the collinear divergences via mass factorization, multiplying by 
T{j. The Z-factors are given by the generic formula: 



Zij(N, a s ,e) 



7ij,o 



:7im,07mj,0 + /?07y,0 



2e 



7ii,i 



g7m,07nm,07mj,0 ' 



A)7ii 



s,07mj,0 + g A)7ij,0 



g (7zm,l7mj,0 



s,07mj,l) + g (A)7y,l 



3e 



which has to be adapted for the various flavour decomposed combinations. 
The indices i,j here either run over i,j £ {l,g} or denote the non-singlet 
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combination. The pure-singlet Z-factor is given by: Z qq = Z qq — Z^s- 
Jij } k are the k-loop anomalous dimensions and denote the expansion 
coefficients of the beta function. The Tij are in the massless case the inverses 
of the Z-factors, however, since we always have at least one heavy quark 
loop, they differ and the Tij contribute at most to two-loop subdiagrams. 
Let us consider, e.g., the renormalized gluonic matrix element Aq 9 , where 
Aq g denotes the mass- and charge-renormalized expression. One finds then: 

= ZgqAQqTqg + Z qq AQ g T g ^ + Z ' ^ Agq^T ^ + Z q g Agg^T 

A (2) 1(2) 7-1,(1)^(1) £-1,(1) Ml) , 7-1,(2) , ( 1(1) , 7 -l, p-l,(l) 

A Qg ~ A Qg + A Q1 A Qg + A gg + ^<?9 + { A Qg + ^qg J 1 gg 

As the general expression in the first line already indicates, there is a mixing 
with Ag q and A g9j Q from 0(af ), while A gq ^Q starts contributing from 0(a 2 s ) 

and therefore at O(af) to AQ g . For Aq 9 given in the second line, one finds 

the 0(e) term of the gluonic one-loop OME, a^, entering the two-loop 
expression via renormalization, as described for example in [3,4,6]. In as 

(3) 

the same way, for the renormalization of Ag g , the 0(e)-term of the two- 
loop expression ag g is needed - as are the terms a gg , aQ q ' PS due to the 

above mentioned operator mixing, and the 0(e 2 ) of the one-loop Aq^. The 
calculation of these 0(e) terms is a first step towards a 3-loop calculation, 
cf. section 3.1., while first calculations of 3-loop diagrams for fixed Mellin 
N are described in section 3.2. 

As a last remark, note that we consider charm quark contributions here, 
while for heavier quarks decoupling [10] has to be applied. 



3.1. Two-loop diagrams to 0(e) for general Mellin N 

Our calculation is performed in Mellin space, where the convolution of 
functions becomes a simple product. The 0(e) terms for the unpolarized 
gluonic OMEs, as for the pure-singlet and non-singlet cases, have been given 
in [6] for general Mellin N. The corresponding polarized contributions are 
to be published soon. The calculation is performed in two ways: on the one 
hand, we rewrote the OMEs in terms of Mellin-Barnes integrals and used M. 
Czakon's mathematica package MB [11] to obtain numeric results, serving 
as a check for the analytic results, which have been obtained expressing the 
OMEs as generalized hypergeometric functions. Expanding these functions 
in e, one has to re-sum the expression for the desired order, which we did 
using integral techniques and C. Schneider's mathematica package SIGMA 
[12]. The results are then given in terms of nested harmonic sums [13, 14], 
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to which we applied algebraic and analytic simplifications [15, 16] to find 
the easiest possible representation. 

(3) 

The term A gg has been newly calculated and is given up to order 0(e) by: 

99, Q 

^{-1 (- 3 l Sl + _,„ 6 « N * + N + v.„) + i (-^ Sl + 



3 3(JV - 1)N(N + 1)(JV + 2) j e\ 9 9{N - 1) AT 2 (N + l) 2 (N + 2) 



8 16(N 2 + N + 1)C2 567V + 47 2P 3 



-C2S1 H — 4 Si 

3 3(JV - 1)JV(JV + 1)(JV + 2) 27{N + 1) 27(N — 1)N 3 (N + l) 3 (N -f- 2) 

■ _8 s _ 20 _ gj 16(JV 2 + TV + lKg P5 

9 3 1 9 3(JV + 1) 9(JV - 1)JV(JV + 1)(JV + 2) 81(N — l)A r4 (A r + 1) 4 (JV + 2) 



4Pi< 2 328 Af 4 + 256 AT 3 - 247 AT 2 - 175 A" + 54 2 AT - 

9(Af - l)Af 2 (Af + l) 2 (Af + 2) 81(JV - 1)JV(JV + l) 2 1 + 3(JV - 

16(AJ 2 + AT + 2) 2 \ if 4P 2 



»} 



+T F C F 



(N - 1)N 2 (N + 1) 2 (N + 2) / £ \ (JV - 1)JV 3 (JV + 1) 3 (JV + 2) 
4(JV 2 + JV + 2) 2 C 2 -P4 



(A" - 1)A" 2 (A" + 1) 2 (A" + 2) (A" - 1)A" 4 (A" + 1) 4 (N + 2) / 

4(AT 2 + AT + 2) 2 C 3 P2C2 Pe 



3(Af - l)AT 2 (Af + 1) 2 (AT + 2) (AT - 1)N 3 (N + 1) 3 (N + 2) 4(JV - 1) JV 5 (JV + l) 5 (JV + 2) 



)} 



Pi = 3A/ a + 9JV 5 + 22JV 4 + 29JV 3 + 41iV 2 + 28JV + 6 , 

P 2 = N s + 4jV 7 + 8iV 6 + 6JV 5 - 3JV 4 - 22iV 3 - 10iV 2 - 8N - 8 , 

P 3 = 15Ar 8 + 60JV 7 + 572JV 6 + 1470Ar 5 + 2135JV 4 + 1794JV 3 + 722JV 2 - 24JV - 72 , 

P 4 = 15iV 10 + 75JV 9 + 112YV 8 + UN 7 - 61JV 6 + 107A/ 5 + 170JV 4 + 36JV 3 - 36iV 2 - 32JV - 16 , 

P 5 = 3A/ 6 + 9iV 5 + 307JV 4 + 599A/ 3 + 746JV 2 + 448iV + 96 , 

P 6 = 1473JV 10 + 7365JV 9 - 40888JV 8 - 201406JV 7 - 377863JV 6 - 391687JV 5 - 230434JV 4 , 

-56896JV 3 + 480JV 2 - 4608JV - 3456 , 

P 7 = 2JV 6 + 4iV 5 + N 4 - ION 3 - 5iV 2 - 4iV - 4 , 

P 8 = 4jV 10 + 16JV 9 + 4jV 8 - 92iV 7 - 156A/ 6 - 74jV 5 + IIjV 4 - 17JV 2 - 12JV-4 . 

where we have calculated an all order e result, which is solely given in terms 
of Euler T and t/j functions. This expression is also needed in the context of 
the variable flavour number scheme []. 

In the unpolarized case, all 2-loop 0(e) terms are now known. In the 
polarized case, the calculation proceeds in the same way and we calculated so 
far the gluonic, pure-singlet and non-singlet terms, which will be published 
soon [6]. 



3.2. Fixed values of N at three loops 

As a next step towards a full 0(af) calculation, we started calculat- 
or 

ing unpolarized three-loop OMEs A^-'q for fixed values of Mellin N. The 
contributing OMEs are: singlet: {A Qg , A g9i Q, A gQi Q}, pure-singlet: Aq^, 
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non-singlet: {A qq q , A qq q , A q}, where we have operator mixing be- 
tween the singlet and pure-singlet terms. The first object of investigation 

(3) 

is the gluonic Aq g : The necessary three-loop diagrams are generated using 
QGRAF [17], where the operator product expansion has been implemented 

up to insertions of operators with three and four gluonic lines. The number 

(3) 

of diagrams contributing to Ag g , e.g., is 1478 diagrams with one and 489 
diagrams with two quark loops, where at least one of the loops is heavy. 
The steps for the calculation of these self-energy type diagrams with one ad- 
ditional scale set by the Mellin variable N, are the following: The diagrams 
are genuinely given as tensor integrals due to the operators contracted with 
the light -cone vector A, A 2 = 0. The idea is, to first undo this contrac- 
tion and to developed a projector, which, applied to the tensor integrals, 
provides the results for the diagrams for a specific (even) Mellin N under 
consideration. So far, we implemented the projector for the first 4 contribut- 
ing Mellin N, N = 2, ...,8, where the color factors are calculated using [18]. 
The diagrams are then translated into a form, which is suitable for the 
program MATAD [19] by M. Steinhauser, which does the expansion in e 
for the remaining massless and massive three-loop tadpole-type diagrams. 
We have implemented these steps into a FORM [20] program and tested it 

(2) 

against some of our two-loop results and the all-order e result of A ' q and 
found agreement. We then turned to a subset of the 3-loop diagrams, the 
diagrams oc Tp: T f Cf, T f Ca, which are currently under investigation. 

4. Conclusions and outlook 

We calculated the 0(e) contributions to heavy flavour Wilson coefficients 
for general Mellin variable N at 0(a 2 ), as a first step towards a O(a^) 
calculation. Furthermore, we installed a program chain to calculate the 
corresponding 3-loop diagrams to O(a^), with the help of MATAD. This 
chain is now existing and we expect first results in the near future. 
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